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ABSTRACT 

The Sloan Lens ACS (SLAGS) Survey is an efficient Hubble Space Telescope (HST) Snapshot imaging 
survey for new galaxy-scale strong gravitational lenses. The targeted lens candidates are selected 
spectroscopically from within the Sloan Digital Sky Survey (SDSS) database of galaxy spectra for 
having multiple nebular emission lines at a redshift significantly higher than that of the SDSS target 
galaxy. The SLAGS survey is optimized to detect bright early-type lens galaxies with faint lensed 
sources, in order to increase the sample of known gravitational lenses suitable for detailed lensing, 
photometric, and dynamical modeling. In this paper, the first in a series on the current results of our 
HST Gycle-13 imaging survey, we present a catalog of 19 newly discovered gravitational lenses, along 
with 9 other observed candidate systems that are either possible lenses, non-lenses, or non-detections. 
The survey efficiency is thus > 68%. We also present Gemini 8m and Magellan 6.5m integral- field 
spectroscopic data for 9 of the SLAGS targets, which further support the lensing interpretation. A 
new method for the effective subtraction of foreground galaxy images to reveal faint background 
features is presented. We show that the SLAGS lens galaxies have colors and ellipticities typical 
of the spectroscopic parent sample from which they are drawn (SDSS luminous red galaxies and 
quiescent main-sample galaxies), but are somewhat brighter and more centrally concentrated. Several 
explanations for the latter bias are suggested. The SLAGS survey provides the first statistically 
significant and homogeneously selected sample of bright early-type lens galaxies, furnishing a powerful 
probe of the structure of early- type galaxies within the half-light radius. The high confirmation rate 
of lenses in the SLAGS survey suggests consideration of spectroscopic lens discovery as an explicit 
science goal of future spectroscopic galaxy surveys. 

Subject headings: gravitational lensing — galaxies: elliptical and lenticular, cD — galaxies: evolution 
— galaxies: formation — galaxies: structure 



1. INTRODUCTION 

In the currently favored cosmological scenario, the 
matter content of the universe is dominated by a 
cold and dark component of unknown particle species 
whose only significant interaction with the smaller bary- 
onic matter fraction occurs through the gravitational 
force. This "cold dark matter" (GDM) picture is most 
strongly required by observations on the largest scales 
HSnergel et a,1.ll2nnf iTegmark et a,1Jl2nn4HPerciva,l et a,1.l 
120021) . The GDM scenario holds that galaxies form 
within the potential wells of extended dark-matter ha- 
los which began their collapse while baryonic matter 
was still ionized ijWhite fc ReesllTgTSt iBlumenthal et all 
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IT981 . This view is supported by direct evidence for 
dark matter in disk galaxies from the observation of ro- 
tational velocities that remain approximately constant 
out to radii at which the stellar g alactic component 
makes a diminishing contrib ution fe.g. lRubin et al.119801: 
lva,n A1ba,da fc Sa,ncisiiri986l). 

Unlike disk galaxies, early-type galaxies (E and SO) are 
generally pressure supported and lack bright kinematic 
tracers at large radii such as Hj. As such their kinemat- 
ics are more difficult to measure and interpret in terms of 
mass density profiles. The density structure of early-type 
galaxies is nevertheless of great interest for numerous rea- 
sons. First, their structure is the "fossil record" o f their 
formation and evolution processes (e.g. Wcchslcr^^D 
12002; Zha o et alJ l2003: Locb & Peebles 2003; Gao et alJ 
12004') . Hierarchical GDM galaxy- formation theories hold 
that early-typ e galaxies are built through the merging 
of late-typ es fKauffmann. White, fc Guiderdonil Il993t 
iBaugh. Cole . & Frenk 1996J, which should have pre- 
dictable consequences for the structure of the merger 
products. The most stringent test of these theories 
will require precise observational measurements of early- 
type mass profiles. For example, detailed measure- 
ment of the structure of high-surface-brightness early- 
type galaxies will enable quantitative tests of the GDM 
theory on scales where baryonic and radiative processes 
have significant effects upon the structure of the host 
dark-matter halo (e.g. through "adiabatic contraction" 
e.g. Florcs ct al. 1993; Blumcnthal et al. 1986; Mo et al] 
119981 iGnedin et al.l l2004l) , altering it significantly rel- 
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ative to the form expected to re sult from pure coUi- 
sionless dark-matt er collapse (e.g. iNavarro et alJ 1199^ 
iMoore et alJll99^ . Second, early- type galaxies exhibit 
great regularity in their photometric, spectroscopic, and 
kinematic properties, as described, e.g., by the well- 
known Fundamental Plane (FP) relation between veloc- 
ity dispersion, effective radius, and surface bri ghtness 
()Diorgovski fc Davis! IT987t Dressier et alJ Il987|) . The 
"tilt" of the FP relative to the simple expectation of 
the virial theorem can be understood in terms of a de- 
pendence of galaxy str ucture and mass-to-light ratio on 
total galaxy ma ss '( Ben der et al.|ll99a iCiotti et al]ll99fil: 
iBertin et alT20n2L.Truiillo et al.l2004D . However, further 
constraints on the mass structure of early-type galaxies 
are needed in order to distinguish between the various 
effects of diff ering stellar populat ions (the explanation 
given by, e.g.. lGerhard et al.l200H) . differing density pro- 
files ("structural nonhomology" or "weak homology"), 
and differing dark-matter fra ctions (the explanation of, 
e.g.. IPadmanabhan et al.|l2004t) in giving rise to the FP. 

Despite great progress, observational results remain 
uncertain because of the small number of individual 
galaxies suitable for study. Stellar-dynamica l measure- 
ments of local earl y-type galaxies (e.g . Bcrtin et al.ll99'^ 
Franx et al.lll994t iGerhard et al.ll2001t iCappellar i et all 
2005ir the statistics of earl y-type gra vitation al-lens 
galax i es (e.g. ^ ochanck 1996; QhSISS; Rusi n et a,].l 
120031 iRusin fc Kochanek .2005^1 . and combined lens- 
ing and dynamical mea surements of the few systems 
amenable to such study (Kooomans & Treu"2002', 1200,'^ 
IXrcu & Kooomans 2002, 2003, 2004, hereafter KT) gen- 
erally argue for the presence of a significant amount of 
dark matter even on the scale of the half-light radius, 
leading to an approximately linear increase of enclosed 
mass with radius and thus to a flat equivalent rotation 
curve as in disk galaxies. However, in some cases the 
presence of dark matter is not required, and the observed 
kinematics can be reproduced with a constant mass-to- 
light ratio (e.g. IBertin et alJ 119941: iRomanowskv et aP 
|2nQ3), corresponding to a declining rotation curve. Fur- 
thermore, [Kochanek (2003) has pointed out an apparent 
conflict between isothermal mass proflles (i.e. flat rota- 
tion curves) and some gravitational-lens time delays un- 
der the assumption of Hq « 70 km s~^ Mpc~^. Due to 
this persistent uncertainty about the structure and diver- 
sity of early- type galaxies (see also Kochanek 2004), it is 
important to collect precision data for a larger number 
of objects, possibly spanning a large range of redshifts, 
environments, and mass. 

Strong gravitational lensing provides the most direct 
probe of mass in early- type galaxies: a measurement 
of the mass enclosed within the Einstein radius. When 
combined with a spatially resolved measurement of the 
line-of-sight velocity-dispersion proflle and the surface 
brightness of the lens galaxy, this mass measurement 
can be used to constrain the luminous and dark-matter 
mass profiles simultaneously through the Jeans equation 
(e.g. KT). Unfortunately, only a handful of known strong 
lenses are amenable to this type of analysis, often because 
the back-ground source is a bright QSO whose images 
outshine the lens galaxy. 

We have therefore initiated the Sloan Lens ACS 
(SLAGS) Survey (see iBolton et aP l2005j) to discover 
a much larger sample of new early-type strong grav- 



itational lenses suitable for detailed photometric and 
dynamical study and thereby realize the full poten- 
tial of gravitational lensing as a probe of early-type 
galaxy structure. S LAGS uses the method described by 
IBolton et alJ (j2004', hereafter B04) to select candidate 
galaxy-scale gravitational lens systems from within the 
Sloan Digital Sky Survey (SDSS) spectroscopic database 
on the basis of multiple higher-redshift emission lines in 
the spectrum of lower-redshift target galaxies. For dis- 
cussion of other spectroscopic lens surveys and lens dis- 
coveries based on the identification of anomalous emis- 
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data also provide stellar velocity dispersion measure- 
ments for all SLAGS lens candidates. The most promis- 
ing candidates (see § I2.1|l are observed through F435W 
and F814W using the Advanced Gamera for Surveys 
(AGS) aboard the Hubble Space Telescope (HST). AGS 
images enable us to confirm bona fide lenses, measure 
detailed photometric and morphological parameters of 
the lensing galaxies, and obtain accurate astrometry and 
surface-brightness measurements of strongly lensed fea- 
tures with which to constrain gravitational-lens models. 
Finally, confirmed lenses from the program are being tar- 
geted for deep ground-based spectroscopy in order to ob- 
tain spatially resolved velocity-dispersion profiles of the 
lensing galaxies, which will provide more precise dynam- 
ical constraints on their mass structure. 

This paper (Paper I) is the first in a series presenting 
the results of the current Gycle-13 HST imaging compo- 
nent of the SLAGS Survey. This paper presents the cat- 
alog of new lenses conflrmed by the SLAGS ffST-AGS 
Snapshot Survey and their properties as measured by 
SDSS, describes our methods of image processing and 
lens-galaxy image subtraction, and provides an analysis 
of the possible selection biases that bear upon whether 
or not our lenses are a representative sample of simi- 
lar SDSS early-type galaxies. We also present ground- 
based integral-fleld spectroscopy of several SLAGS sys- 
tems which supports a strong-lensing interpretation of 
the observed features. Paper II presents photometric and 
morphological measurements of the SLAGS lens sample 
from (HST) imaging and places the sample within the 
context of the FP, Paper HI presents gravitational-lens 
and dynamical modeling results, and Paper IV focuses 
on the properties of the lensed source galaxies. 

The outline of this paper is as follows: In §13 the survey 
selection procedure and the available data of each sys- 
tem is described. In §|21 nineteen newly discovered lens 
systems and nine unconfirmed systems or non-detections 
are presented. In §0] the statistics and possible selec- 
tion effects of the sample are discussed. §|5l discusses 
some of the implications of the SLAGS survey for fu- 
ture surveys and in §0 conclusions are drawn and future 
work is discussed. Throughout this paper, we assume 
a cosmological model with flM = 0.3, J^a = 0.7, and 
Ho = 70/i7okms"i Mpc^^ (with /170 = 1). 

2. THE SURVEY 

2.1. Candidate Selection 

The method by which we select our lens candidates 
is essentially that described by B04. Briefly, we sub- 
tract best-flt PGA templates (a byproduct of the redshift 
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TABLE 1 

Properties of Observed Systems 



System Name 



Plate-MJD-Fiber 



RcS (") 



2FG -^BG fa (kms ^) Sample Lens 



SDSS 
SDSS 
SDSS 
SDSS 
SDSS 
SDSS 
SDSS 
SDSS 
SDSS 
SDSS 
SDSS 
SDSS 
SDSS 
SDSS 
SDSS 
SDSS 
SDSS 
SDSS 
SDSS 
SDSS 
SDSS 
SDSS 
SDSS 
SDSS 
SDSS 
SDSS 
SDSS 
SDSS 



J003753. 
J021652. 
J073728. 
J081931. 
J091205. 
J095320. 
J095629. 
J095944. 
J102551. 
J111739. 
J120540. 
J125028. 
J125135. 
J125919. 
J133045. 
J140228. 
J142015. 
J154731. 
J161843. 
J162746. 
J163028. 
J163602. 
J170216. 
J171837. 
J230053. 
J230321. 
J232120. 
J234728. 



92^ 
30^ 

42^ 



21-094220.1 
54-081345.3 
45+321618.5 
453444.8 
002901.1 
520543.7 
77+510006.6 
07+041017.0 
31-003517.4 
60+053413.9 
43+491029.3 
25+052349.0 
70-020805.1 
05+613408.6 
53-014841.6 



-632133.5 
-601914.8 
-572000.0 
-435327.4 
005357.5 
-452036.2 
-470729.5 
-332044.7 
-642452.2 
-002237.9 
-142217.9 
093910.2 
000521.2 



0655- 
0668- 
0541- 
0441- 
0472- 
0902- 
0902- 
0572- 
0272- 
0835- 
0969- 
0847- 
0337- 
0783- 
0910- 
0605- 
0788- 
0617- 
0815- 
0364- 
0626- 
0627- 
0973- 
0352- 
0677- 
0743- 
0645- 
0684- 



52162-392 
52162-428 
51959-145 
51868-108 
51955-429 
52409-577 
52409-068 
52289-495 
51941-151 
52326-571 
52442-134 
52426-549 
51997-480 
52325-279 
52377-503 
52353-503 
52338-605 
52072-561 
52374-337 
52000-084 
52057-518 
52144-464 
52426-464 
51789-563 
52606-520 
52262-304 
52203-517 
52523-311 



18.00. 
19.07. 
19.38, 
18.63. 
17.31, 
18.57, 
18.41. 
18.52, 
17.07. 
18.72, 
18.40. 
18.40. 
18.58. 
18.80, 
18.34. 
18.26, 
16.39. 
17.94, 
18.78. 
18.54. 
18.84. 
18.92. 
18.01, 
16.86. 
18.97. 
17.58. 
16.07. 
19.81. 



16.81 
17.46 
17.84, 
17.51 
16.22 
17.61 
17.17, 
17.48 
16.03 
17.56 
17.22 
17.26 
17.59 
17.46 
17.45 
16.98 
15.56 
16. 
17.60 
17.29 
17.41 
17.68 
16.91 
15.97, 
17.63 
16.39 
15.21 
18.51 



,16.39 
,16.90 
,17.15 
,17.07 
,15.78 
,17.22 
,16.62 
,17.02 
,15.57 
,17.12 
,16.65 
,16.77 
,17.24 
,17.01 
,17.05 
,16.49 
,15.17 
,16.39 
,17.09 
,16.89 
,16.92 
,17.18 
,16.40 
,15.54 
,17.14 
,15.96 
,14.82 
,17.95 



2.16 ± 0.06 
3.05 ± 0.13 

2.16 ± 0.13 

2.32 ± 0.13 
3.36 ± 0.05 

1.77 + 0.09 

2.33 ± 0.09 
1.21 ± 0.04 
4.05 ± 0.08 
2.49 ± 0.11 
2.30 + 0.10 
1.76 ± 0.07 
3.64 + 0.19 
1.94 + 0.07 
0.84 + 0.04 
2.67 + 0.08 

2.17 + 0.03 
2.56 + 0.06 

1.34 + 0.05 
2.08 + 0.08 
2.02 + 0.07 
1.48 + 0.05 
2.80 + 0.07 
3.67 + 0.07 
1.76 + 0.10 
3.02 + 0.09 
3.92 + 0.05 

1.78 + 0.28 



0.1954 
0.3317 
0.3223 
0.1943 
0.1642 
0.1310 
0.2405 
0.1260 
0.1589 
0.2285 
0.2150 
0.2318 
0.2243 
0.2333 
0.0808 
0.2046 
0.0629 
0.1883 
0.1989 
0.2076 
0.2479 
0.2282 
0.1784 
0.0899 
0.2285 
0.1553 
0.0819 
0.4168 



0.6322 
0.5235 
0.5812 
0.4462 
0.3239 
0.4670 
0.4700 
0.5350 
0.2764 
0.8232 
0.4808 
0.7946 
0.7843 
0.4488 
0.7115 
0.4814 
0.5350 
0.3955 
0.6656 
0.5241 
0.7933 
0.6745 
0.4357 
0.7367 
0.4635 
0.5170 
0.5324 
0.7145 



265 + 10 
332 + 23 
310 + 15 
231 + 16 
313 + 12 
207 + 14 
299 + 16 
212 + 12 
247 + 11 
279 + 21 

235 + 10 
254 + 14 
216 + 23 
263 + 17 
178 + 09 
275 + 15 
194 + 05 
243 + 11 
257 + 25 
275 + 12 
260 + 16 
221 + 15 
239 + 14 
270 + 16 
283 + 18 
260 + 15 

236 + 07 
330 + 50 



LRG 

LRG 

LRG 
MAIN 

LRG 
MAIN 

LRG 
MAIN 

LRG 
MAIN 
MAIN 
MAIN 
MAIN 

LRG 
MAIN 

LRG 
MAIN 

LRG 
MAIN 

LRG 

LRG 
MAIN 

LRG 
MAIN 

LRG 

LRG 
MAIN 

LRG 



Yes 
Yes 
Yes 
? 

Yes 
? 

Yes 
Yes 

? 

7 
Yes 
Yes 
Yes 

7 

Yes 
Yes 
Yes 

7 

Yes 
Yes 
Yes 

7 

7 
Yes 
Yes 
Yes 
Yes 
7 



Note. — System Name gives truncated J2000 RA and Doc in the format HHMMSS.ss+DDMMSS.s. De Vaucoulcurs model SDSS (AB) mag- 
nitudes have been de-reddened using dust maps from lScti lcgel. Finkbcincr. & Davis (1998), and have statistical errors of appproximatcly 0.01. De 
Vaucouleurs effective radii are quoted at the intermediate axis. Parent Sample column indicates whether candidate was selected from the SDSS 
luminous red galaxy sample or from the SDSS main galaxy sample (sec ^ 12. 11 . 



pipeline) from the observed SDSS target-galaxy spec- 
tra. We require the SDSS continuum to be well-fit by 
the template, which effectively yields a parent sample of 
galaxies with well-behaved absorption-dominated spec- 
tra and very secure redshifts, which we denote zpQ (for 
"foreground"). The residual spectra are then scanned 
for nebular line emission at redshifts greater than zpc- 
Spectra in which such emission is significantly detected 
for at least three separate common atomic transitions 
at a single background redshift (zbg) are taken as lens 
candidates for having emission at two different redshifts 
along the same line of sight as sampled by the 3-arcsec 
diameter SDSS spectroscopic fiber. 

Taking zpG; -zbGj and the measured stellar velocity dis- 
persion (Ta from SDSS allows us to determine to first or- 
der the strong-lensing cross section for each system in 
the source plane - using a singular isothermal sphere 
(SIS) model - as ttO"^, with the Einstein radius given by 
9e = 4:'K{a^/c'^){DLs/Ds)- In the image plane the re- 
gion of multiple imaging is four times larger than the 
source-plane strong-lensing cross section, in the sense 
that all images within a radius 26 e of the potential cen- 
ter will show strongly-lensed counterimages. To max- 
imize the number of strong lenses in our survey, we 
formed our Cycle-13 HST target list from the candidates 
with the highest predicted lensing cross section. Thus 
our HST target sample is velocity-dispersion selected to 
leading order, with an additional selection bias in favor 
of systems with larger angular-diameter distance ratios 
Dls/Ds (a function of the redshifts zbg and zfg)- We 
take 20 of our targets from the list of lens candidates 



published in B04, which were selected from within the 
luminous red galaxy spectros copic sample of the SDSS 
(LRG. lEisenstein et alJl200l|) . The LRG sample is de- 
fined by photometric selection cuts that very efficiently 
select massive early-type galaxies in the redshift range 
0.15 < zpG ^ 0.5, as confirmed by SDSS spectroscopy. 
These galaxies are very homogeneous in their spectral, 
photometric, and morphological properties, and we place 
no further requirements on our LRG-sample lens can- 
di dates beyond the i r pho tometric selection as described 
in lEisenstein et all lj2001 f) and subsequent spectroscopic 
confirmation as galaxies with redshifts z > 0.15. The re- 
maining 29 targets on our Snapshot list are selected with 
the same spectroscopic algori thm from within the MAIN 
galaxy sample of the SDSS ijStrauss et al.ll2002() . The 
MAIN sample is much more heterogeneous, and we im- 
pose a quiescent, absorption-dominated spectral criterion 
by requiring our lens candidates to have rest-frame equiv- 
alent widths in Ha of EWna <1.5A (with a few excep- 
tions made in the interest of maximizing the total num- 
ber discovered lenses in the program: see SDSSJ1251- 
0208 below.). Although most of the systems we target 
do indeed exhibit early-type morphology, the selection of 
SLACS lens candidates from within the LRG and MAIN 
galaxy samples of the SDSS is purely spectroscopic. 

2.2. ACS Image Processing 

Our ACS Snapsh ot observing strategy is discussed in 
IBolton ct"ar (2005) , and consists of one 420s Wide-Field 
Channel (WFC) exposure through each of the two fil- 
ters F435W and F814W. The SLACS lens candidates 
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Fig. 1. — SDSS spectroscopy showing background line emission in SLACS target galaxies. Upper gray line shows SDSS template fitted 
to the continuum of the foreground galaxy, and lower gray lines shows the l-cr noise level. 



are selected to yield bright lenses with faint background 
sources. This facilitates the photometric, morphological, 
and kinematic study of the lens galaxy, but can also make 
the relatively faint lensed features difficult to detect and 
even more difficult to use in lens modeling. The key to 
success is effective subtraction of the image of the lens 
galaxy by fitting a smooth model to the image data of 
the c andidate lens galaxy (e.g. Pens et al. 2002: Simard 
Il99i^ . Given the extremely regular isophotal structure 
of most early-type galaxies, this is a reasonable proposi- 
tion. The most common parameterized model for early- 
type brightness distributions is the generalized de Vau- 



couleurs or Sersic law l|Sersicll968HCiottll99lh . We first 
attempted to use elliptical Sersic model-fitting and sub- 
traction in our analysis, but generally found that the 
global systematic residuals of the fit were large com- 
pared to the surface-brightness levels of the lensed fea- 
tures that we hope to use to constrain gravitational-lens 
models. Significant galaxy-core residuals are especially 
pronounced in the Sersic-subtracted F814W data. The 
shortcomings of the Sersic model led us to implement a 
more generalized galaxy-model fitting procedure involv- 
ing a b-spline fit to the radial profile with a low-order 
multipole dependence to fit the angular structure, which 
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we describe in Appendix IXI 

3. NEW LENSES AND OTHER OBSERVED SYSTEMS 

Here we present the catalog of 28 candidate lens sys- 
tems^ observed by the SLACS survey through 2005 
March 22. Photometric and spectroscopic parameters 
for the sample as measured by SDSS are given in Ta- 
ble n Figure ^ presents the SDSS discovery spectra of 
all targeted systems, focused on the background-redshift 
line emission. For definiteness, we enumerate the several 
possible explanations for any one of our spectroscopically 
selected lens candidates as follows: (1) a multiply imaged 
background galaxy ("strong lens"), (2) a singly imaged 
though possibly magnified background galaxy (a "non- 
lens" in our classification), (3) a multiple-image system 
corresponding to a projection of multiple singly imaged 
sources with similar colors, and (4) spurious noise fea- 
tures in the SDSS spectrum. Given the adopted signif- 
icance threshold for selection and the careful treatment 
of the noise properties of SDSS spectra (B04), we ex- 
pect (4) to be an unlikely explanation, though it mer- 
its consideration in the case of any lower signal-to-noise 
candidates that do not produce detections in follow-up 
observations. Systems falling into category (3) could in 
principle be mistaken for bona fide strong lenses, but 
can be expected to have ill-fitting lens models, since real 
lenses occupy only a small subspace of all conceivable im- 
age configurations. For simplicity, this paper only cat- 
egorizes ACS-observed systems into "strong lenses" on 
one hand and "everything else" on the other, with "ev- 
erything else" including systems for which a definitive 
explanation cannot be made based on the data in hand 
(i.e. ACS non-detections). 

To determine the incidence of strong lensing, we ex- 
amine the F814W and F435W residual images formed 
by subtracting smooth models of the foreground galaxy 
constructed as described in Appendix El If the resid- 
ual image of a target shows multiple images with sim- 
ilar colors that can be reproduced by a simple lens 
model with the potential center fixed at the optical cen- 
ter of the foreground galaxy we classify the system as 
a strong gravitational lens. The lens models for these 
systems are presented in Paper III, but we summarize 
the lens- modeling process here. The pr ocedure is based 
on the method of lWarre n &: Dy j i2003tl. with an im ple- 
mentation a s desc ribed in lTreu fc Koopman^ l)2004|) and 
iKoonman^ 1)20051) . We parametri ze the lens galaxy as a 
singular isothermal ellipsoid fSIE. lKormann et al.l ll994\ 
with three free parameters corresponding to lens strength 
(i.e. Einstein radius), ellipticity, and major-axis position 
angle. The source-plane brightness distribution of the 
candidate lensed galaxy is represented on a pixel grid. 
We determine the best-fit lens parameters and source- 
plane pixel brightnesses by minimizing ^ with respect 
to the galaxy-subtracted /-band image-plane data, while 
also imposing a brightness-dependent regularization on 
the source-plane distribution to suppress curvature. Un- 
successful strong-lens models are characterized by an in- 
ability to map putative multiply-imaged features to a sin- 
gle region of the source plane for any set of SIE model pa- 
rameters. SDSSJ1025 is an example of a system with an 

* The Cycle— 13 program is still ongoing and another 11 systems 
have been observed, with 9 more being scheduled. 



apparent arc and counterimage which nonetheless can- 
not be reproduced by an SIE lens model, and is thus not 
classified as a strong lens. ACS data for target galaxies 
classified as lenses are shown in Figure |21 We note again 
that all target systems have line emission in their SDSS 
spectra at a redshift higher than that of the foreground 
galaxy, and thus the lensing interpretation does not rest 
on imaging data alone. 

Figure 01 presents all systems other than those that we 
classify solidly as lenses. This figure encompasses sys- 
tems that may indeed be strong lenses as well as systems 
that appear to be definite Konlenses. We note that a 
significant number of systems that we consider possible 
but not definite lenses show a faint candidate counterim- 
age near the center of the foreground galaxy opposite a 
more prominent image at larger radius. If these features 
were all due to simple residual error in the foreground 
galaxy subtraction, this configuration would not be ex- 
pected. The features could conceivably be explained 
by the prominent image "pulling" the foreground-galaxy 
model to the side, leading to an under-subtraction of 
the foreground-galaxy fiux near the center on the oppo- 
site side. We believe this explanation is unlikely since 
we have taken care to mask neighboring images when 
fitting the foreground-galaxy model (Appendix^), and 
we do not see any corresponding over-subtraction on 
the near side of the core. The most asymmetric dou- 
ble lenses will in general have faint counterimages at 
small angular offset from the lens center and will be the 
most difficult lenses to confirm, although their status as 
lenses or non-lenses can have a significant impact on sta- 
tistical inferences based on the lens sample. Integral- 
field spectroscopy provides the best chance to associate 
background-redshift line emission with these faint pu- 
tative counterimages near the lens center and thereby 
confirm or reject a strong-lensing hypothesis. 

Many of the lens galaxies shown in Figure [3 show such 
striking features in their residual images that their in- 
terpretation as strong lenses is effectively certain when 
considered in combination with their SDSS spectroscopic 
detection and gravitational- lens models (Paper III). For 
putative lenses with less dramatic morphology, further 
evidence as to their status may be obtained with integral- 
field spectroscopy, which can confirm or deny the spatial 
coincidence between the high-redshift line emission and 
the candidate lensed features seen with HST. For several 
of our SLACS target galaxies, we have obtained spatially 
resolved spectroscopy with the integral-field units (IFUs) 
of the Ina mori Magellan Area! Cam era and Spectrograph 
(IMACS, IRigelowfcnress1eHl2003l ISchmol] et a1Jl200l 
on the 6.5m Walter Baade (Magellan I) telescope at 
Las Campanas Observatory and the 8m Gemini- North 
Multi-Obje ct Spectrograph fCMOS-N. fHook et al.i,20 03: 
iMiirrav et al. 2003) at Mauna Kea; we present these data 
in Figure El In all cases, the IFU emission-line images 
are coincident with the ACS residual-image morphology, 
supporting the strong-lensing interpretation. 

We note here that the background sources in our lens 
systems are generally faint galaxies with irregular mor- 
phology. As evidenced by the emission-line flux by which 
they are originally detected, they should also be regarded 
as starforming. The absence of significant broadening 
of the Balmer emission lines at the ^ 150km s~^ res- 
olution of the SDSS further suggests that they do not 
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host quasars. For a recent observational study of the 
likely source population of emission- line galaxies, see 
iDrozdovskv et all l|2005|) . 

4. STATISTICS AND SELECTION EFFECTS 

The statistics of strong gravitational lensing are 
the subject of a great body of litera ture (e.g. 
iTurner et al."1984'^ 'Fu kugita fc Turned IT99ll iKochanekl 
^96; Chac ct al. 2002), most of it dealing with the lens- 
ing of quasars by galaxies. The rigorous statistics of 
the SLACS survey, though involving similar considera- 
tions, would be sufficiently distinct to require their own 
detailed treatment: our original sample is formed from 
potential lensing objects rather than from potentially 
lensed sources, and we must consider the effects of an 
extended source image and finite fiber sampling. We de- 
fer a full lens-statistical analysis of the SLACS survey to 
a future paper. Here we address the selection effects of 
the SLACS survey in terms of the robust statistics of the 
distribution of lens-galaxy observables. 

As discussed above, the SLACS target selection in- 
volves an explicit velocity-dispersion and redshift selec- 
tion. This should not in itself present any complication in 
the comparison of our lenses to other galaxies with sim- 
ilar velocity dispersions and redshifts. After controlling 
for this selection, we may ask whether our lenses are a 
representative sample of the populations from which they 
have been drawn: luminous red galaxies and quiescent 
main-sample galaxies from the SDSS. If so, their relation 
to early-type galaxies selected by other means is simply 
determined by the SPSS LRG a nd MAIN sa mple target 
select ion function ijEisenstein et al. 2001; St rauss et al.l 
120021) . The question can be framed between two al- 
ternatives: do we select simply for spectral superposi- 
tion caused by a chance foreground-background coinci- 
dence, or is there some significant bias by which galaxies 
with certain properties are more likely to be selected as 
lens candidates based on their SDSS spectra? To test 
this possibility we exploit the parent sample of ^ 10^ 
SDSS galaxies (Early Data Release through Data Re- 
lease 3). For each lens we construct a control sam- 
ple of SDSS galaxies with the same redshift and ve- 
locity dispersion (within the uncertainties) — quantities 
for which we explicitly select — and test whether the 
lens galaxy is typical of galaxies in the control sample 
in terms of its magnitude, color, effective radius, and 
isophotal ellipticity — quantities for which we do not ex- 
plicitly select. Although this is not a definitive test of 
the representative nature of our lenses (it will not un- 
cover any "hidden variable" bias), it is straightforward, 
robust, and informative. By forming control samples 
at the redshift of the lenses, we avoid the necessity of 
applying evolution- and /c-corrections to the broadband 
magnitudes. We may perform our comparison sensibly 
for galaxies with velocity dispersions well-measured by 
SDSS. We exclude the lens SDSSJ1251-0208 from the 
analysis, both because its spectral signal-to-noise ratio is 
too low for a confident velocity-dispersion measurement 
and because its Ha equivalent width exceeds our quies- 
cent threshold. (In fact, it is a bulge-dominated spiral- 
galaxy lens). The double-lens galaxies SDSSJ1618+4353 
and SDSSJ1718+6424 are also excluded, as is the lens 
SDSSJ1205-t-4910, which has significant flux from a 
neighboring galaxy within the SDSS spectroscopic fiber. 



This leaves a well-defined sample of 15 single early-type 
lenses. The comparison sample for each lens consists of 
unique galaxies from the SDSS database with redshifts 
within ±0.005 of the lens redshift, velocity dispersions 
within ±15kms^^ (approximately the median velocity- 
dispersion error of the sample) of the lens galaxy, and 
median signal-to-noise per spectral pixel greater than 
8. For LRG-sample targets, we also require the com- 
parison sample to pass the photometric LRG cuts, and 
for the MAIN-sample lenses, we require a rest-frame 
EWna < 1-5 A. We make a luminosity-distance cor- 
rection to the broadband magnitudes and an angular- 
diameter-distance correction to the effective radii within 
the redshift slice to place all galaxies closer to the exact 
redshift of the lens. The resulting individual comparison 
samples have as few as 17 and as many as 1793 galaxies, 
with a mean of 451 and a median of 337 control galaxies 
per lens galaxy. 

We would like to answer the question of whether our 
lenses have observables (luminosities, colors, effective 
radii, and ellipticities) consistent with having been drawn 
from the distribution of those observables seen in the 
control samples. Since we only have one lens at each 
redshift, a straightforward Kolmogorov-Smirnov (K-S) 
test of the lens sample against the control samples is 
not possible: even under the null hypothesis, each indi- 
vidual lens has a different parent distribution. We may 
however put all these parent distributions on common 
footing by using the fact that the K-S test is invariant 
under a monotonic rescaling of the variable under con- 
sideration. Specifically, for each individual lens's control 
sample, there exists a transformation of the observable of 
interest that converts the control distribution into a uni- 
form distribution between the minimum and maximum 
values. The correspondingly transformed lens observable 
is simply equal to the normalized rank (between and 
1) of the lens quantity within the cumulative distribu- 
tion of the control sample in that quantity. Thus we may 
perform a one-dimensional K-S test of the distribution of 
lens-observable ranks against a uniform distribution over 
the interval to 1 . This test is in some sense like a rank- 
correlation test within the K-S formalism: does ranking 
as a lens correlate with ranking in luminosity, color, size, 
or ellipticity? 

Figure [51 shows the cumulative distributions of lens 
rankings in r (magnitude) , g ~r (color) , Re (effective ra- 
dius), b/a (isophotal axis ratio), and le (r-band effective 
surface brightness: le oc 10~^''^^ R^^) within their con- 
trol samples. From these distributions we compute the 
statistic -Dks j equal to the maximum difference between 
the cumulative rank distribution and the null-hypothesis 
distribution. Since the parent distribution under the null 
hypothesis is known exactly by construction (uniform 
probability for any rank between and 1), and since the 
total number of control galaxies is much larger than the 
number of lenses, we compute the significance of Dks 
for a distributio n of 15 values against a known parent 
distribution (e.g. iPress^t al...l992'l . (The statistical sig- 
nificance of an outlying lens will be limited by a smaller 
control sample in that its rank is at least 1/iVgai and 
at most 1 — 1/A^gai, where -A^gai is the number of galax- 
ies in its control sample.) The resulting probabilities of 
the lens sample having been drawn at random from the 
control-sample distributions in r, g — r, i?e, b/a, and 
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Fig. 2. — HST ACS-WFC F814W and F435W imaging of new gravitational lenses from the SLACS survey. Also shown are residual 
images with smooth b-spline lens-galaxy models subtracted, revealing lensed features more clearly. Images are formed from flat-fielded 
single-image native ACS data. Cosmic-ray and other zero-weight pixels are replaced with values from a median-smoothed residual image, 
with the b-spline model galaxy values added for the direct images. Images are 8" X 8", with N up and E left. All systems in this figure are 
modeled successfully with singular isothermal ellipsoid lens models, as described in Paper III. 



are 0.097, 0.655, 0.264, 0.550, and 0.085 respectively. In 
terms of color and axis ratio, the SLACS lenses seem 
to be a representative sample of the parent distribution 
of SDSS galaxies. However, the distribution of SLACS 
galaxies in brightness and effective radius (and hence in 
surface brightness) is somewhat significantly skewed to- 
wards brighter and more centrally concentrated systems. 

This trend is in the same sense as that of the FP: at 
fixed velocity dispersion, brighter galaxies have smaller 
effective radii. Paper II shows that the SLACS lenses 
fall on the local FP corrected for stellar evolution, and 
further examines the degree of lens bias towards higher 
effective surface brightness in the context of the FP. Al- 
though the statistical significance is not great, the trend 
is at first glance suggestive of the effect that more cen- 
trally cond ensed objects are m ore efficient gravitational 
lenses (e.g. iLi fc Ostrikeill2002() . to the extent that stel- 
lar mass constitutes the dominant lensing component. A 



larger sample of similarly selected lenses will allow us to 
test this selection effect with greater statistical signifi- 
cance. Possible interpretations of this finding — as dis- 
cussed in detail in Papers II and III — include observa- 
tional selection effects due to the finite size of the SDSS 
fiber which may bias our survey towards the highest sur- 
face brightness lens galaxies. In Paper III we will see 
that, within the context of powerlaw mass models, the 
SLACS lens sample shows great homogeneity in the slope 
of the total mass density profile (7') and no significant 
correlation between 7' and Jo. This uniform degree of 
mass concetration within the lens sample itself argues 
somewhat against a lens-selection bias due to mass con- 
centration, but this argument is limited by the lack of 
an equally direct probe of the mass distribution in the 
non-lens control sample. 

5. DISCUSSION AND IMPLICATIONS FOR FUTURE 
SURVEYS 
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Fig. 2. — (continued) 

The large number of new strong lenses confirmed by 
the SLACS Survey demonstrates the technical feasibility 
of carrying out spectroscopic surveys for strong gravita- 
tional lenses, and suggests that similar gravitational lens 
surveys should be considered as explicit goals of future 
large-scale spectroscopic surveys. The SLACS results arc 
particularly notable when one considers that this type of 
lens survey was never one of the science or design goals of 
the SDSS, and the survey is thus a testament to the per- 
formance of the SDSS telescope and spectrograph, and 



to the optimal quality of the SDSS spectroscopic data 
pipeline. With the clarity of hindsight, we may identify 
those particular features of the SDSS that were benefi- 
cial to the spectroscopic lens survey project and those 
that could conceivably be improved in any future survey 
to increase the yield of spectroscopically selected strong 
lenses suitable for specified science goals. 

The large parent sample of galaxies observed is a clear 
advantage afforded by the SDSS for this project or any 
other rare object survey. Empirically, we have found that 
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(continued) 



spectroscopic lenses and lens candidates appear with a 
frequency on the order of one in a thousand in the SDSS: 
starting with ^ 10^ early-type galaxies is therefore one 
key to success. At the same time, the large number of 
highly magnified galaxies that we have found close to 
the spectroscopic detection limit (note the generally low 
emission-line signal-to-noise ratio in the discovery spec- 
tra of Figure Q suggests that we may only be seeing 
the tip of the iceberg. We conclude that (1) a spectro- 
scopic lens survey will always be at an advantage if it 



can "piggy back" on a more broadly conceived spectro- 
scopic survey that does not have lenses as its sole scien- 
tific goal, and (2) future spectroscopic surveys that go 
significantly deeper than the SDSS should, all other fac- 
tors being equal, discover a significantly larger fraction 
of gravitational lenses. 

One might suspect that a significant factor for SLACS 
survey is the large (3-arcsec diameter) aperture of the 
SDSS spectroscopic fibers. Certainly a larger fiber will 
always collect more photons, but for lenses with image 
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Fig. 2. — (continued) 

splittings of ~ 2" (i.e. ~ 1" Einstein radii) such as we 
have presented here, a 3-arcsec fiber could afford an ad- 
vantage over smaller fibers disproportionate to its larger 
aperture duo to the concentration of Icnscd flux away 
from the center of the fiber. We can test the importance 
of this aperture-size effect by simulating observations 
with varying fiber apertures. We assume an Einstein 
radius of 1'.'2 (an approximate median value from Paper 
III) and an exponential-disk background galaxy with a 
disk scale length of 0'.'2. We compute lensed images for 



a range of source-galaxy impact parameters using a sin- 
gular isothermal sphere lens model. We then smear the 
images with 1'.'3 seeing (the median SDSS spectroscopic 

seeing for our targets) and integrate over the 3-arcscc 
diameter SDSS fiber. If we reduce the fiber aperture 
to 2-arcscc in diameter, we find that the flux collected 
from the lensed galaxy is reduced to a fraction of 0.4-0.5, 
depending on impact parameter. This is essentially the 
factor of (1/1.5)^ ~ 0.44 corresponding to the diminished 
aperture size. Thus for lenses of the angular size typical 
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Fig. 3. — HST ACS-WFC imaging of lens candidates observed by the SLACS survey for which the ACS data arc inconclusive as to the 
incidence of strong lensing. Images are as in Figure |^ See notes on individual systems in Appendix El 



of SLAGS lenses, the advantage of a large fiber is just the 
simple advantage of a large aperture due to the smearing 
effects of seeing; a smaller fiber would not have lost a 
disproportionate amount of line flux due to lensing ge- 
ometry. However, for lenses with larger image splittings 
(which would occur for higher redshift sources behind 
the same foreground galaxies), a small fiber would likely 
be a qualitative disadvantage to a spectroscopic lens sur- 
vey. Because a fiber-diameter much greater than 3 arcsec 
would probably increase spurious non-lensed interlopers, 
we conclude that the SDSS fiber-size represents a very 
good compromise, when searching for massive early-type 
lens galaxies at these redshifts. 

For two principal reasons, our spectroscopic survey 
owes much success to having targeted lens candidates 
with background redshifts confirmed by multiple emis- 
sion lines. First, the incidence of false-positive spectro- 
scopic detections and emission-line mis-identifications is 
negligible. Second, to fixed line-flux limits, intermediate- 
redshift emission from [O ii] 3727, [O ill] 5007, and 



the hydrogen Balmer series is more common than 
Lyman-g emission at high redshift ijHippelein et al. 20031: 
iMaier et alJl200l . and thus we see an abundance of 
oxygen/Balmer lenses. However, for optical surveys, 
this survey strategy limits the background redshifts to 
zbg S 0.8, beyond which H/? and [O ill] 5007 move 
out of the observable band. There is a definite incen- 
tive to discover significant numbers of lenses with higher 
lens and source redshifts, both to probe evolution in the 
lens population and to observe lensed images at larger 
physical radii within the lens galaxy in order to obtain 
greater leverage on dark-matter halos. Therefore in de- 
signing future surveys to discover gravitationally lensed 
emission-line galaxies at higher redshift, two obvious con- 
siderations are increased survey depth to detect a fainter 
population, and increased spectroscopic resolution be- 
yond the A/ A A « 1800 resolution of the SDSS in order 
to split the [O ii] 3727 doublet and resolve the charac- 
teristic asymmetric profile of Lyman-a emission (thus 
permitting more secure single-line redshifts). 
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Fig. 3. — (continued) 



One of the greatest uncertainties at the outset of the 
SLAGS survey was simply the fraction of strong lenses 
within the candidate list. One can attempt to calcu- 
late the lensing probabilities of the candidates given zfGj 
zbg, CTa, assuming a model for the SDSS spectroscopic 
observations, and assuming a luminosity function (LF) 
and size/shape for the background em ission-line galax- 
ies, as was done in lBolton et al.l l)2004j) . but the ingredi- 
ents are all uncertain. To put it simply, we did not know 
beforehand whether or not the majority of our targets 



would owe their high-redshift line emission to the PSF- 
smeared wings of bright unlensed galaxies at large impact 
parameter rather than to Einstein-ring images of faint 
galaxies closely aligned with the SDSS target galaxy. The 
high ratios of lenses to non-lenses and of quads/rings to 
doubles that we see are suggestive of a large magnifi- 
cation bias in our selection, with highly magnified faint 
lensed galaxies detected with greater frequency than less 
magnified lensed galaxies or unlensed projections. Future 
work will quantify the extent of any such magnification 
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Fig. 4. — IFU narrowband imaging of several SLACS target systems. Left panels show HST-ACS F814W residual images, for reference. 
Center panels show IFU emission-line images, and right panels show IFU continuum images. IFU images are constructed from IFU spectra 
by fitting linear continuum and Gaussian emission-line models as described in Bolton &: Buries (in preparation). Data are from IMACS 
(SDSSJ0037, SDSSJ0216, and SDSSJ2321) and CMOS-N (others). Emission-line images for SDSSJ0216, SDSSJ0956, and SDSSJ1547, have 
been smoothed spatially with a 7-lenslet hexagonal kernel to suppress noise. Narrowband images are formed at the redshifted wavelength of 
the following background-galaxy emission Unes: [O in] 5007 (0037, 0737, 1402, 1702), (0956), and [O ii] 3727 (0216, 1547, 1630, 2321) 

bias, and will derive statistical constraints on the lens 
and source populations based on the incidence of lensing 
within the survey. 

6. CONCLUSIONS AND FUTURE WORK 

The HST Cycle-13 SLACS Survey^ has produced a 
catalog of 19 previously unknown early- type strong grav- 
itational lens galaxies. The > 68% fraction of genuine 
lenses within the the first 28 candidates observed shows 
that the survey strategy — spectroscopic candidate selec- 
tion from within the SDSS, followed by HST-ACS Snap- 
shot observations of the systems with the largest esti- 
mated lensing cross sections — is an efficient means for the 
discovery of new gravitational lenses; similar lens surveys 
should be considered as an explicit goal of future spec- 
troscopic surveys. By targeting galaxies with the high- 
est estimated lensing cross sections, the Cycle-13 SLACS 
survey has effectively selected high-mass lensing galaxies 
(strong lensing cross section is simply proportional to 

8 HSTprogram #10174; PI: Koopmans. 




lens mass within the Einstein radius). The natural selec- 
tion of lenses in other surveys is also weighted by lensing 
cross section, and thus strong lenses are typically mas- 
sive galaxies. During HST Cycle-14, we will extend our 
survey to galaxies of lower mass in order to use strong 
lensing and stellar dynamics to measure the mass depen- 
dence of early-type galaxy structure and mass-to-light 
ratio^". These galaxies will have a lower lensing cross 
section and hence should have a lower lensing rate, al- 
though magnification bias may skew the distribution of 
candidates in favor of lenses over non-lenses. The final 
combined lens sample will be a unique resource for the 
detailed measurement of the mass profile of early type 
galaxies within the effective radius. 

The SLACS sample represents the largest single cat- 
alog of uniformly selected early-type gravitational lens 
galaxies assembled to date. The lensed images pro- 
vide an aperture mass constraint within a typical scale 
of ~ i?B/2. The SLACS lens galaxies are all brighter 

1° HSTprogram #10587; PI: Bolton; 118 Snapshot targets. 
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Fig. 4. — (continued) 

than their lensed background galaxies by many magni- 
tudes, and are thus ideally suited to detailed photomet- 
ric measurement as with the current ACS data. They 
are also excellent candidates for spatially resolved spec- 
troscopy to constrain dynamical galaxy models in combi- 
nation with the mass constraints from strong lensing. In 
fact, nearly all of the SLACS galaxies already have well- 
measured luminosity- weighted velocity dispersions inside 
the 3-arcsec SDSS fiber (see Table and thus we have 
already nearly tripled the number of known gravitational 
lenses with known stellar velocity dispersions. Paper III 
will present the constraints on the logarithmic slope of 
the total radial density profile (luminous plus dark mat- 
ter) that can be obtained by combining lens models with 
the aperture-integrated SDSS velocity dispersions. We 
may derive tighter constraints and a significant decompo- 
sition into luminous and dark components by obtaining 
spatially resolved line-of-sight velocity dispersion mea- 
surements within the lens galaxies (see KT). 

An analysis of the distribution of SLACS lens-galaxy 
photometric and structural parameters within control 
samples from the SDSS database shows that the SLACS 
lenses are typical of their SDSS parent sample with re- 
gard to color and ellipticity. The location of SLACS lens 
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galaxies within the FP is the subject of Paper II. How- 
ever, the SLACS lens galaxies exhibit a somewhat signif- 
icant bias towards brighter and higher surface-brightness 
galaxies at fixed redshift and velocity dispersion (§0) see 
also paper II). The significance of this bias can be con- 
firmed with a larger sample of similarly selected lenses, 
which we plan to obtain with the continuation of our 
ACS survey during EST Cycle-14. 

Results in this paper are based in part on data from 
the Sloan Digital Sky Survey (SDSS) archive. Funding 
for the creation and distribution of the SDSS Archive 
has been provided by the Alfred P. Sloan Foundation, 
the Participating Institutions, the National Aeronautics 
and Space Administration, the National Science Foun- 
dation, the U.S. Department of Energy, the Japanese 
Monbukagakusho, and the Max Planck Society. The 
SDSS Web site is http://www.sdss.org/ The SDSS 
is managed by the Astrophysical Research Consortium 
(ARC) for the Participating Institutions. The Participat- 
ing Institutions are The University of Chicago, Fermilab, 
the Institute for Advanced Study, the Japan Participa- 
tion Group, The Johns Hopkins University, Los Alamos 
National Laboratory, the Max-Planck-Institute for As- 
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Fig. 4. — (continued) 

tronomy (MPIA), the Max-Planck-Institute for Astro- 
physics (MPA) , New Mexico State University, University 
of Pittsburgh, Princeton University, the United States 
Naval Observatory, and the University of Washington. 

These results are also based in part on observations 
obtained with the 6.5-ni Walter Baade telescope of the 
Magellan Consortium at Las Campanas Observatory. 

These results are also based in part on observations 
obtained under program GN-2004A-Q-5 at the Gemini 
Observatory, which is operated by the Association of 
Universities for Research in Astronomy, Inc., under a 
cooperative agreement with the NSF on behalf of the 
Gemini partnership: the National Science Foundation 
(United States) , the Particle Physics and Astronomy Re- 
search Council (United Kingdom) , the National Research 
Council (Canada), CONICYT (Chile), the Austrahan 
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Research Council (Australia), CNPq (Brazil) and CON- 
ICET (Argentina). 

These results are also based in part on observations 
made with the NASA/ESA Hubble Space Telescope, ob- 
tained at the Space Telescope Science Institute (STScI), 
which is operated by the Association of Universities for 
Research in Astronomy, Inc., under NASA contract NAS 
5-26555. These observations are associated with program 
#10174. Support for program #10174 was provided by 
NASA through a grant from STScI. 

The work of LAM was carried out at Jet Propulsion 
Laboratory, California Institute of Technology, under a 
contract with NASA. 

The authors thank the anonymous referee for a careful 
reading of the manuscript and for a constructive report. 



APPENDIX 

A. B-SPLINE GALAXY MODEL SUBTRACTION 

The b-spline technique is a well-known method for fitting a piecewise-defincd polynomial of arbitrary order to the 
dependence of a series of data values upon an independent variable (e.g. dc Boor 1977). The coefficients of the 
polynomial change at breakpoints in the independent-variable domain, whose spacing may be chosen to allow more or 
less freedom depending upon the level of detail to be fit. A b-spline of order n (where by convention rt = 1 is piecewise 
constant, n = 2 is linear, n = 3 is quadratic, n = 4 is cubic, and so on) has continuous derivatives to order n — 2 across 
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Fig. 5. — Kolmogorov-Smirnov tests of the rank of lens-galaxy observables within the cumulative distributions of those observables 
in control galaxy samples with similar redshifts and velocity dispersions. From upper left to lower right are the distributions for r-band 
magnitude, g — r color, effective (half-light) radius Re, isophotal axis ratio b/a, and effective surface brightness /e (all as measured from 
SDSS imaging data for a de Vaucouleurs galaxy model). Solid black lines show the cumulative distribution in rank, solid gray lines show the 
null hypothesis for this distribution, and vertical dash-dot lines show the location of the maximum difference between the two distributions. 

the breakpoints. For our galaxy modeling we use b-splines of order n = 4, which have continuous 0th, 1st, and 2nd 
derivatives. The coefficients of the b-spline are determined by fitting to the data in a least-squares sense, given the 
breakpoint spacing and derivative-continuity conditions. A b-spline model may be reinterpreted in terms of a number 
of localized basis functions within the domain whose shapes are set by the order of the b-spline and by the breakpoint 
spacing and whose amplitudes arc determined by the fit to the data. The fit itself is entirely linear, and the localized 
infiuence of each basis function within the domain implies that only a banded-diagonal matrix need be inverted in the 
solution for the coefficients. 

The radial b-spline technique for galaxy images permits smooth fitting of arbitrary radial brightness profiles. Since 
we do not have multiple dithered exposures to combine, we perform our b-spline galaxy-model fitting in the native 
ACS pixel coordinates of the images. We account for the distortion in the ACS by using the solution provided in 
the image headers to compute relative tangent-plane RA and Dec values for all pixels in the image, which are taken 
as the independent variables for the fit. We fit for surface brightness as a function of position using cosmic-ray- 
masked flat-fielded images (i.e. calibrated by the CALACS software pipeline^^ to measure surface brightness rather 
than counts-per-pixel) , so our fits are not biased by slight variations in pixel area across the images. The one non- 
linear step in our b-spline galaxy-model fitting is the determination of the center of the lensing galaxy; we describe 
our centering method further below. The model-fitting procedure is carried out using the adopted center for the lens 
galaxy within a suitable subsection of the ACS field (typically 12" x 12"), along with an error image, a mask specifying 
cosmic-ray and other zero-weight pixels, and a mask corresponding to stars, neighboring galaxies, and any apparent 
backg round-galaxy features. The cosmic-ray masks are initially generated by the LACOSMIC software (|Sn Dokkui3 
1200 1|) and adjusted manually over the image subsection surrounding the lens; the neighboring-object masks are created 
manually. For each pixel in the sub-image, the radial offset R from the galaxy center (in arcseconds) and the azimuthal 
angle 9 relative to a fixed position angle are computed. A set of breakpoints in R is chosen for the fit, typically every 
0'.'2-0'.'3 in the central 1" to 2", with increased spacing further out. (We note again that the ACS pixel scale is 0'.'05 
per pixel.) For the case of complete Einstein rings, breakpoints at the ring radius can be removed to allow smoother 
interpolation of the lens-galaxy model over the ring region. A multipole angular dependence is incorporated into the 
fit as follows. If the purely one-dimensional b-spline fit is represented as 

/(i?)=^afc/fc(i?) , (Al) 
where the fk (R) are the localized basis functions that arc only non-zero over a small range of breakpoints and the 

(see|http : //g««. stsci . edu/hst/acs/ j 
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are their amplitudes, then the two-dimensional fit is represented by 

I{R, 0)=Y, [b,nk cos{m9) + Cmk sin(m0)] fk (R) . (A2) 

The number of multipole orders to be fit is chosen individually for each galaxy, with m — (monopole) and m — 2 
(quadrupole) always present and higher-order terms added if necessitated by systematic angular structure in the 
residual images. With the occasional exception of an m = 1 (dipole) term to capture slightly disturbed morphologies 
or mild mis-centering in the central regions, only even multipole orders are used. This form of angular dependence 
is ideal for fitting early-type galaxies: the fit remains linear, the global symmetries seen in early-type galaxies are 
naturally reflected in the low-order terms, and effects such as isophotal twists, varying ellipticity with radius, and 
diskiness/boxyness can be captured with minimal effort (with only the last of these effects requiring multipole orders 
beyond the quadrupole). We note also that the monopole term can automatically include a fit to the sky background. 

Following the completion of the initial b-splinc model- fitting step (which for our 12" x 12", 240 x 240-pixel images 
generally takes 2-3 seconds on a 2.53GHz Pentium 4 Linux PC), we examine the residual (data— model) image for 
faint features not associated with the smooth galaxy model and not initially masked, and perform a second fit with an 
updated mask, also adding multipole orders as needed. Since at this stage we are principally concerned with generating 
high-quality residual images and not with the measurement of lens-galaxy parameters, we fit directly to the images 
without convolving the b-spline model with the HST point-spread function. As an example, we describe the fit to 
the sky-subtracted image of the E3 lens galaxy SDSSJ0912+0029 to give the F814 residual image shown in Figure |21 
which uses breakpoints spaced every 0"3 in R and multipole terms of order 0, 2, and 4 (monopole, quadrupole, and 

octopole). The ratio of the m — 2 amplitudes (defined as ^^2*; + '^2fe) to = amplitudes rises from ^ 0.25 in 

the center to « 0.4 at the effective radius (about 3"), and the ratio of the to = 4 to the to = amplitudes rises from 
« to « 0.1 over the same range. Some of this increase in angular structure with radius is due to our fitting directly 
to the PSF-blurred data. 

The centering necessary for the b-spline galaxy fit is accomplished by first fitting with respect to a best-guess center 
with monopole, dipole, and quadrupole terms. This fit is then evaluated in a ring and used to determine a fiux centroid, 
which is adopted as a new best-guess center for another iteration of the same procedure. Within a few iterations this 
process — which amounts to minimizing the dipole term — converges to a position that we adopt as the constrained 
center for the final model fit. 

B. NOTES ON INDIVIDUAL SYSTEMS 
B.l. Lenses 

SPSS J0216— 0813: A faint but definite counterimage is seen to the NE in both bands in this system, opposite the more 
obvious extended arc. 

SDSSJ0912+0029: This system shows a faint extended counter-arc to the S in the F814W band, opposite the more 
prominent Northern arc. Both images are only marginally detected in the FASbW band. 

SDSSJ0956+5100: This system includes a small round red companion to the N of the lens galaxy, not to be confused 
with the prominent Icnsed images. 

SDSSJ1251— 0208: This galaxy has extended spiral structure in addition to the lensing bulge. 

SDSSJ1330— 0148: This system is judged to be a lens on the basis of the compact counterimage detected at high 
significance in both bands. 

SDSSJ1402-I-6321: This lens , with faint quadruple images, was the first system observed by the SLACS survey and is 
the subject of lBolton etai\ l|2005D . 

SDSSJ1618-I-4353: The angular resolution of HST reveals this system to be a pair of foreground galaxies. The pair 
lenses a compact background source into a 3-f 1 quad configuration. The extra compact image to the west of the 
counterimage is more blue in color than the lensed images and thus does not pertain to the lensed configuration. 

SDSSJ1718-I-6424: A second nearby galaxy contributes significantly to the lensing potential in this system. 

B.2. Other Systems 

SDSSJ1117+0534, SDSSJ12594'6134: The imaging depth is not great enough to unambiguously show whether the 
putative lensed features appear in both bands, but the general geometry is plausibly similar to J0912-I-0029, which 
has an extremely solid gravitational lens model (Paper III). 

SDSSJ1636-f4707: The curvature of the arcs is suggestive of ring geometry, and the knot-counterknot features along 
the North-South axis may be images of each other. Deeper observations would reveal whether or not this system is a 
real Einstein ring, but it is not certain from the current data. 

SDSSJ1702-I-3320: Both F435W and F814W data show possible counter-image features with plausible relative orien- 
tations and separations at low signal-to-noise. 
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